Magnetic order in the Shastry-Sutherland model 
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The ground state properties of the Shastry-Sutherland model in the presence of an external field 
are investigated by means of variational states built up from unpaired spins (monomers) and singlet 
pairs of spins (dimers). The minimum of the energy is characterized by specific monomer-dimer 
configurations, which visualize the magnetic order in the sectors with fixed magnetization M = 
S/N. A change in the magnetic order is observed if the frustrating coupling a exceeds a critical 
value a c (M), which depends on M. Special attention is paid to the ground state configurations at 
M = 1/4, 1/6, 1/8. 
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I. INTRODUCTION 



The Shastry-Sutherland modefl defined by the two- 
dimensional spin 1/2 Hamiltonian 



H 



E 

<x,y> 



S(x)S(y) + a ^ S(x)S(y) (1.1) 

«x,y» 



with nearest neighbor couplings and frustrating next- 
nearest neighbor couplings on the diagonals shown in Fig. 
|l| has attracted a lot of interest for theoretical and exper- 
imental reasons: 




FIG. 1. The couplings in the Shastry-Sutherland model. 
Nearest and next-nearest neighbor couplings are represented 
by dotted and solid lines. 

(1) The product wave function 

*= n ^ 

((x,y» 

built up from singlet states 

[x,y] = i=(x + (x) x _(y)-x-(x)x + (y)) (1.3) 



is known to be an eigenstate of the Hamiltonian (1.1), 
which turns out to be the ground state, if the coupling 



a exceeds a critical value a c (a c w 1.4).0 The phase dia- 
gram has been studied recently by Weihong, Oitmaa and 
HamerB by means of series expansions. 

(2) The Hamiltonian (1.1) is suggested to be an appro- 
priate model for the compound SrCu2(BOa)2 the mag- 
netic properties of which have beexLjinvestigated in recent 
high magnetic field experiments.aQ Plateaus have been 
found in the magnetization curve M = M [B) at ratio- 
nal values of the magnetization MjMg = 1/3, 1/4, 1/8, 
where Mg — 1/2 is the saturating magnetization. 

From the theoretical point of view the appearance 
of magnetization plateaus is well understood in quasi 
one-dimensional systems, e.g. with ladder geometry.Q 
Here, a quantization rule has been formulated by Os- 
hikawa, Yamanaka and Afft£ck,tj which originates from 
the prediction of soft mpdesEI based on the Lieb-Schultz- 
Mattis (LSM) theorem]!!] Only the position of the pos- 
sible plateaus - i.e. the quantized value of the magneti- 
zation - are predicted by this rule. The upper and lower 
critical fields which define the width of the plateau, how- 
ever, depend on the magnitude of transition matrix el- 
ements with a momentum transfer corresponding to the 
relevant soft mode. These matrix elements contribute to 
the dynamical and static structure factors and a strong 
peak in these quantities at the soft mode momenta is 
needed for a pronounced plateau in the magnetization 
curve. El 

The extension of the Lieb Schultz Mattis construction 
to higher dimensions (D >-X) meets difficulties. As was 
pointed out by OshikawaHl magnetization plateaus are 
possible in higher dimensions as well, provided that the 
"commensurability condition" is satisfied. Based on a 
topological argument he shows, that this condition is a 
robust non-pcrturbative constraint. 

The emergence of magnetization plateaus in a modified 
Shastry-Sutherland model has been discussed in Ref. [ 
Recently, Misguich, Joliccer, and Girvin studied the 
emergence of magnetization plateaus in the framework of 
a Chern-Simons theoryt^l 

In Ref. [ |J] Kageyama etal. proposed that product 
wave functions of the type fll.2| ) with certain distributions 
of Ns singlets (1.3) and Nt triplets 
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{x,y} = x+(x)x+(y) 



(1.4) 



might yield an appropriate ansatz for the ground state 
in the sector with total spin 



S T = N T /A 



(1.5) 



where Nt and N$ are constrained by the total number 
of sites 



N S + N T = N/2 . 



(1.6) 



Typical examples of these states are shown in Figs. ||(b), 
||(b), |](c) and^(d). An effective Hamiltonian, describing 
the interaction between singlets and triplets, has been 
developed and evaluated in Refs. [ p^ , p!5[ . 

In this paper, we investigate a wider class of prod- 
uct wave functions - which we call monomer-dimcr 
configurations!!! - and which arc aimed to describe 
the magnetic order at those magnetizations (M — 
1/4, 1/6, 1/8, 1/16), where plateaus are expected. Indeed 
we find a change in the magnetic order if the frustration 
parameter exceeds a critical value a c (M), which depends 
on the magnetization M. For a > a c (M) we recover the 
singlet-triplet configurations proposed in Ref.[ |4|]. For 
a < a c (M), however, we find new configurations with 
lower energy. 

The outline of the paper is as follows: In Sec. II we 
define the monomer-dimcr configurations. In Sec. Ill 
we m inimize the expectation value of the Hamiltonian 
(1.1) between monomer-dimer configurations. This pro- 
cedure singles out specific configurations, which visual- 
ize the magnetic order at fixed magnetization M. In 
Sec. IV we introduce the "frozen monomer approxima- 



tion" , which allows to lower the energy expectation values 
between monomer-dimer configurations without chang- 
ing the magnetic order, i.e. the distribution of "frozen" 
monomers. 

The quality of the frozen monomer approximation is 
studied for M = 1/8, 1/6, 1/4 in Sec. |V| by a compari- 
son with the ground state energies obtained from exact 
diagonalizations on finite clusters. 

We also look for the formation of magnetization 
plateaus. Possible interpretations of the observed 
plateaus in SrCu2(BO$)i are discussed in Sec. VI. 



II. MONOMER-DIMER CONFIGURATIONS 



The Shastry-Sutherland Hamiltonian (1.1) conserves 
the total spin 



and we therefore start from eigenstates of 



5(5 + 1) and 5 3 



-5,. 



,5. 



(2.1) 



(2.2) 



Following HulthenO, these states \K, v = 25) can be 
constructed in the sector with total spin 5 -i.e. magne- 
tization M = S/N- as product states of 



• unpaired spin- up states at sites xi, . . . , x„, v = 25: 
|x+)=X+(x) (2.3) 
which we call "monomers" 



• singlets of paired spins [x,y] ( |l.3|) at sites x, y 
("dimers"): 

u=2S 

\K, V )= n K+) n^-y] ( 2 - 4 ) 

j =1 <x.y> 

Note, that in the monomer-dimer configuration K each 
site x is occupied exactly once: either by a monomer or a 
dimer. Moreover, monomer-dimer configurations \K,v) 
yield an overcomplete non-orthogonal set of eigenstates 
with total spin 5. 



The expectation value of the Hamiltonian (1.1) be- 
tween monomer-dimer configurations can be easily cal- 
culated with the following rules: 



1 

4 



(2.5) 

(2.6) 
(2.7) 
(2.8) 



(l+,2+|S(l)S(2)|l+,2+) 

([1,2]|S(1)S(2)|[1,2]} = - 
([1,3]2+ |S(1)S(2)|[1,3]2+) = 
([1,3][2,4]|S(1)S(2)|[1,3][2,4]> = 0. 
If we count on each configuration the numbers 

• n[ 0) (K) of nearest neighbor dimers 

• N2°\k) of next-nearest neighbor dimers (corre- 
sponding to Fig. |l|) 

• N^\k) of nearest neighbor monomer pairs 

• (K) of next-nearest neighbor monomer pairs 
(corresponding to Fig. [I]) 

we can immediately compute the expectation value: 
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{K,v\H\K^) = --A° ) -l aN 2 ) 



3 

■ — i 

4 

4 1 4 2 



(2.9) 



In order to minimize this expectation value we have to 
look for configurations with 

• a maximum number of nearest neighbor dimers 
N^ 0) if a < 1 

• a maximum number of next-nearest neighbor 
dimers if a > 1 

• a minimum number of monomer pairs 

on nearest and next-nearest neighbor sites. 
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III. MAGNETIC ORDERING AT FIXED 
MAGNETIZATIONS 

A. M — 1/4, v = N/2 



changes its sign for a — 4/5, which means there is 
a change in the magnetic order from configuration K\ 
to K 2 if the frustration parameter exceeds the value 
a c (M = 1/4) = 4/5. 



Let us start with M — 1/4. In this situation we 
have to distribute v — N/2 monomers and N/4 dimcrs 
on the square lattice. We cannot avoid the appearance 
of monomer pairs on nearest and next-nearest neighbor 
sites, but we can minimize their numbers n[ and N 2 
if we cover the lattice in the way shown in Fig.||(a). The 
bold lines symbolize the nearest-neighbor singlets, the 
thin lines the monomer pairs on nearest and next-nearest 
neighbor sites. Dimer pairs can interact via the next- 
nearest neighbor couplings in the Hamilton ian: they are 
represented by dotted lines. According to (2.£), the ex- 
pectation value of the Hamiltonian is found to be: 



£(#1,1/4) = (Jf 1 ,JV/2|lT|A'i,JV/2) 
N, 

= — (-l + a/4). 



(3.1) 



A second configuration, shown in Fig. y(b) has been 
proposed in Ref . [ p| as a possible ground state configura- 
tion for M = 1/4. In this case, the expectation value of 
H turns out to be: 



E(K 2 ,l/4) = (K 2 ,N/2\H\K 2 ,N/2) 
N 




(b) 

























X 



FIG. 2. Monomer-dimer configurations with minimal en- 
ergy expectation value at M = 1/4. Bold lines represent 
"dimers" i.e. paired spins coupled to singlets. Dotted lines 
indic ate t he couplings between dimers induced by Hamilto- 
nian (1.1). Unpaired spin-up states ("monomers") [Eq.(2.3)] 



are symbolized by solid points; their couplings on nearest and 
next-nearest neighbor sites are indicated by thin lines, (a) 
configuration K\ for a < a c (M — 1/4) = 4/5; (b) configura- 
tion Ki for a > a c (M = 1/4) 



The difference of (3.1) and (p 



N 5 

E{K U 1/4) - E(K 2 , 1/4) = -(-1 + -a) 



B. M = 1/6, v = N/3 

Next, we turn to the case M = 1/6, where we have 
to distribute v = N/3 monomers and iV/3 dimers on the 
lattice. The configuration K\ [Fig. ||(a)] minimizes the 

number N^ °f monomer pairs on next-nearest neigh- 
bor sites, whereas in the configuration K 2 [Fig. ||(b)] 
the next-nearest neighbor sites of Fig. [j] are occupied 
with singlets and triplets in the spirit of Ref.[ |5|. The 
difference of the expectation values of H 



£(#i,l/6) = {K u N/Z\H\K u N/?>) 

4 ^ 12 ; 
E(K 2 , 1/6) = (K 2 ,N/3\H\K 2 , N/3) 
N5a 
~ "IT 



N 11 

E{K X , 1/6) - E(K 2 , 1/6) = - T (l - -a) 



(3.4) 
(3.5) 

(3.6) 



changes sign for a c (M = 1/6) = 12/11. Again we ob- 
serve a change in the magnetic order from Ki to K 2 if a 
passes this value. 



(3.2) 
/ 

/ ■ K • t : 



(a) 



(b) 

/ . / \ v / . / 

\ , 



\ 
\ 



/ 
/ 



(3.3) 



FIG. 3. Same as Fig| 
Ki a < a c (M = 1/6) 
a > a c (M = 1/6) 



for M = 1/6: (a) configuration 
= 12/11; (b) configuration K2 for 
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It is remarkable to note, that in both cases a < 
a c (M = 1/6) and a > a c (M = 1/6) a stripe order of 
the monomers is predicted. 



C. M= 1/8, v = N/4 

In the case of M — 1/8 we have to distribute v — N / 4 
monomers and 37V/ 8 dimers on the square lattice. We 
can avoid now completely the appearance of monomer 
pairs on nearest and next-nearest neighbor sites as is 
demonstrated by the configuration K\ shown in Fig. 
§(a). Owing to the stripe structure, we can also construct 
a second configuration [Fig. f|(b)] with — N/8 

dimers on nearest neighbor sites and = N/ 4 dimers 
on next-nearest neighbor sites. Two further configura- 
tions K 3 and K± have been proposed in Refs. [ p|.[l5|, 



which only contain n[ 0) = 3N/8 dimers and iV{ iJ = N/8 
monomer pairs on next-nearest neighbor sites. 
The corresponding energy expectation values are 



(i) 



E{K x ,l/8) = --N 
E(K 2) l/8) = -^N(l + 2a) 

N 

E(K 3 ,l/8) = E(K 4 ,l/8) = -—a. 



(3.7) 
(3.8) 
(3.9) 



Comparing the expectation values (^)-(|3j|) we expect 
a change in the magnetic order with a: 

a<l: E(Kt) < E(K 2 ) < E{K 3 ) = E(K 4 ) (3.10) 
K a < I : E(K 2 ) < E(Kx) < E(K 3 ) = E(K 4 ) (3.11) 

I < a : E(K 3 ) = E(K 4 ) < E{K 2 ) < E{K X ) . (3.12) 



IV. THE FROZEN MONOMER 
APPROXIMATION (FMA) 

The monomer-dimer configurations which we devel- 
oped in the last section to describe the magnetic order 
in the Shastry-Sutherland model are no t ei genstates of 
the Hamiltonian (1.1). Application of (1.1) onto these 



states will generate new states. In this section we study 
the impact of those couplings in the Hamiltonian, which 
generate interactions only between dimer pairs, i.e. we 
consider an approximation where the v monomers are 
frozen at sites xi . . .x„ in the configuration \K, v). For 
each of these configurations, we define a decomposition 
of the Hamiltonian in three parts: 



H = H V + H(K) + H(xi,K) 



(4.1) 



(a) 



(b) 



• N 



• N 



XT ' 



(c) 



X H X X X X 



\ 



X X X X X X 

///'//'/ 



(d) 

\ , \ j-N , \ 



V. Zt K /XV. 

x 



X 



x , \ 



\ ' X 



V 
/ 



where 



FIG. 4. Same as Fig^ for M = 1/8: (a) configuration 
K\, a < 1; (b) configuration K2, 1 < a < 3/2; (c) and (d) 
configurations Kz and K4, a > 3/2 



a) H u contains all the nearest and next-nearest neigh- 
bor couplings between the sites xi . . . x^, where the 
monomers are located. All other sites are occupied 
with dimers. They form an antiferromagnetic clus- 
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ter K, which are represented by the dimers and the 
dotted connections between dimers in Figs. §(a)- 
@(d). 

b) The cluster Hamiltonians H(K) is defined by the 
nearest and next-nearest neighbor couplings on the 
dimer cluster K. 

c) The nearest and next-nearest neighbor couplings 
in H(pii,K) take into account the remaining inter- 
actions between the monomer at site and the 
dimers in the cluster K. 

The ground state energy E(K) of the antiferromag- 
netic cluster Hamiltonian H(K) 



H(K)f(K) = E{K)^(K) . 



(4.2) 



is obviously lower than the expectation value of H(K) 
between the dimer product wave function on the cluster 
K. 

The product ansatz 



v=2S 



(4.3) 



yields an eigenfunction of H v + H(K) with energy 

E(K, v) = \n1\K) + jN?\K) + E(K) (4.4) 

which again represents an upper bound for the exact 
ground state energy Eq(M — v/2N) of (4.1) in the sector 
with magnetization M = S/N: 



E {M = v/2N) <E{K,v). 



(4.5) 



In the derivation of (4.5) one has to use the fact, that 
the expectation value of t he i nteraction term H(u, K) 
between the product state (4.3) vanishes, since 



(*{K)\Si(y)\*{K)) = yeK, i = 1,2,3. (4.6) 



V. NUMERICAL RESULTS 

In order to check the quality of the frozen monomer 
appr oximation (FMA), we have computed the energies 
(4.4) and compared with exact diagonalizations of the 
Shastry-Sutherland Hamiltonian at fixed magnetization 
M = v/2N on lattices with N = 4 x 4 = 16 and N = 
4 x 6 = 24. 

The strongest effects due to the frozen monomer ap- 
proximation occur at small magnetizations. We therefore 
start with M = 1/8. 



A. M = 1/8 

From Fig. @(a)-|(c) we see that the interactions be- 
tween the dimers generate a two-dimensional cluster 
which contains all dimers in the configuration. In con- 
trast, the dimers in Fig. ^(d) form quasi-one-dimensional 
"stripe" clusters. 



?5 
n 

s 

II 




a 

FIG. 5. Ground state energies per site Eo(N, M, a)/N for 
Shastry-Sutherland lattices of7V = 4x4, 4x6 sites and corre- 
sponding FMA energies of configurations Ki, ...,Kn at mag- 
netization M = 1/8. 

In Fig. H the expectation values E(Kj,M — 1/8, a), 
j = 1, 2, 3, 4 in the frozen monomer approximation - cor- 
responding to the configurations Fig. [l](a)-[|(d) - are rep- 
resented by dotted lines. 

The following points should be noted: 

• The transition in the magnetic order from config- 
uration K2 [Fig. f|(b)] to configuration K3 [Fig. 
f|(c)] occurs here at a larger value of the frustration 
parameter 



a c (M = 1/8) ~ 2.3. 



(5.1) 



At this value the difference in the expectation val- 
ues E{K X ,M = 1/8, a) - E(K 3 ,M = 1/8, a) 
changes sign. For smaller values of a the distri- 
bution of monomers according to Fig. |^(a),(b) 
is favored in comparison with the distribution of 
triplets in Fig. ||(c),(d). 

For a > 1.3 the expectation values 



E(K 3 ,M) = E(K A ,M) 



(5.2) 



corresponding to configurations K3 and K4 [Fig. 
|](c) and Fig. Q(d)] coincide in the frozen monomer 
approxi mat ion. Indeed, here, the dimer product 
ansatz (2.4) is an eigenstat e o f the antiferromag- 
netic cluster Hamiltonian fl4.2Q . In other words: 
The interactions between the dimers [dotted lines 
in Fig. [|(c),(d)] do not lower the ground state ex- 
pectation value. 
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• The expectation values E(Kj,M = 1/8) deviate 
significantly for a < 1.3 from the exact results given 
by the solid curves. Therefore, other distributions 
of monomers should play an important role in the 
exact ground state. 

• For small a, the eact results show a linear behavior 
which is well reproduced in a perturbative expan- 
sion in a: 



iV- 1 S (Af,a) = ei(Af) + — e 2 (M) 



where 



ei(l/8) 
£2(1/8) 



0|S(x)S(y)|0) j 

-0.59 
-0.43 



1,2 



(5.3) 



(5.4) 



are the ground state expectation values of the near- 
est neighbor (j — 1, (x, y) ) and next-nearest neigh- 
bor (j — 2, ((x, y))) spin-spin correlators of the 
unfrustrated Hamiltonian H(a = 0) = Hi(l, 1) 

• Finite-size effects are small, as can be seen from a 
comparison of the exact results for the two systems 
N = 16 and TV = 24. 



-0.2 



a-~a- l *.~n— o-..a.. 




a 



FIG. 6. Ground state energy per site Eo(N, M, a)/N for a 
Shastry-Sutherland lattice of N = 4 x 6 sites and correspond- 
ing FMA energies of configurations Ki , K2 at magnetization 
M = 1/6. 



4 
3 



[1,2] [3,4] -> $(1,2,3,4) 



(5.6) 



B. M = 1/6 

In this case the interactions between the dimers form 
quasi-one-dimensional clusters with stripe geometry as 
can be seen from Fig. |](a),(b). The expectation values 
E(Kj,M = 1/6, a) j = 1,2 in the frozen monomer ap- 
proximation are shown in Fig. ^. The transition point in 
the magnetic order is found her at 



a c (M = 1/6) = 1.2. 



(5.5) 



At this point the exact result of the N — 4 x 6 = 24 
system (solid line) has its maximum. The expectation 
values E(Kj,M = 1/6) deviate significantly for a < 1.2 
from the exact results given by the solid curve. 



C. M = 1/4 

The configuration K\ in Fig. ^(a) is built up from 
4-point singlet clusters. In the frozen monomer approx- 
imation we lower the energy if we substitute each dimcr 
pair 



1 

FIG. 7. The 2-dimer cluster in the configuration K\ 
[Figfa)] 



by the ground state of the 4-point cluster computed 
from the 4-point Hamiltonian 

H(l, 2, 3, 4) = S(1)S(2) + S(3)S(4) + aS(2)S(3) . (5.7) 
The corresponding ground state energy 

E(l, 2, 3, 4) = - _ I (4 - 2a + o?f /2 < -\ (5.8) 

is lower than the energy of the dimer pair. Taking into 
account this effect in the expectation value (3.1) we find 



E(K X ,M= 1/4) = — (4-2a + a 2 
16 



1/2 



(5.9) 



Note that there are no interactions between the dimers in 
the configuration K2 [Fig. |^(b)]. Therefore the ground 
state energy (^^) cannot be lowered through the fr ozen 
mon ome r approximation. The energy differences of ( |5.9| ) 
and (3.2) changes its sign at 



a c (M = 1/4) = -(-1 + v^3) ~ 0.869 . (5.10) 

Below this value the expectation value E{K\,M = 1/4) 
is a very poor approximation for the exact ground state 
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E (N=4x6,a,M=l/4)/N 



2.5 



a 



FIG. 8. Ground state energies per site Eo(N,M,a)/N for 
Shastry-Sutherland lattices ofAf = 4x4,4x6 sites and cor- 
responding FMA energies of configurations K\ , , K2 at mag- 
netization M = 1/4. 



energy, indicating that the true ground state is not repro- 
duced adequately by the frozen monomer approximation. 

In Fig. ^| we com pare the energy expectation values 
E(K U M = 1/4, a) (U) and, E(K 2 ,M = 1/4, a) 
with the exact diagonalization on finite systems with N = 
4 x 4 = 16 and N = 4 x 6 = 24 sites. The maximum of 
the exact ground state energy Eq(M = 1/4, a) is found 
at a ~ 1.5 far beyond the transition point (5.10) from 
configuration K\ to K 2 . 



We have also studied the formation of plateaus in the 
magnetization curve at M = 1/8, 1/6, 1/4 by exact diag- 
onalizations on the finite clusters with N — 4 x 4 = 16 
and N = 4 x 6 = 24 sites. The lower and upper critical 
fields 

B L {M,a) — E a (M, a) — E a (M — 2/N,a) (5.11) 
Bu{M,a) = E (M + 2/N,a) - E (M, a) (5.12) 

were computed from ground state energies Eq(M — 
2/N, a), E (M, a), E {M + 2/N, a) with neighboring to- 
tal spins 5-1,5,5+1, S = M ■ N. The results are 
shown in Fig. | (a) for M = 1/8 (b) for M = 1/6 (c) for 
M = 1/4. 

For a < 1.2 all critical fields are rather a- independent. 
The finite-size effects indicate that the plateau width 

A(M,a) = Bu(M,a) - B L (M,a) for a < 1.2 (5.13) 

will vanish in the thermodynamic limit N — » 00 as it is 
known for the unfrustrated model (a = 0). 

For a — 1.5 the lower critical fields B^iM = 1/8, a) 
and Bl(M = 1/6, a) have a pronounced minimum; be- 
yond this value (a > 1.5) all lower and upper critical 
fields for M = 1/8, 1/6 increase with a. 

For M = 1/8 [Fig. [91(a)] and a > 1.5, finite-size effects 
appear to be small for Bl(M = 1/8, a) but large for 
Bjj(M = 1/8, a). We suggest that the rectangular ge- 
ometry of the 4x6 system might be responsible for this 
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0.5 



1.5 




a 



FIG. 9. Upper(Bij) and lower(_B_L) critical fields for mag- 
netizations M = 1/8 (a), 1/6 (b), 1/4 (c) calculated on Shas- 
try-Sutherland lattices of N = 4 x 4 (a,c) and N = 4 x 6 (a-c) 
sites. 



failure. It breaks the rotational invariance and therefore 
does not allow for the rotated patterns in Fig. [|. 

For M = 1/4 [Fig. ^(c)] and a > 1.5 we observe a 
rather clean signal for the opening of a magnetization 
plateau. 
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VI. DISCUSSION AND PERSPECTIVES 

In this paper, we have investigated the magnetic order 
of the Shastry-Sutherland model at fixed magnetizations 
M = 1/8, 1/6, 1/4. For large enough values of the frus- 
tration parameter 

a > a c (M = 1/8) = 2.3, a > a c {M = 1/6) = 1.2, 
a > a c (M = 1/4) = 0.89 

configurations built up from singlets and triplets on the 
Shastry-Sutherland lattice [cf. FigsJ|(b), |(b), |(c),(d)] 
- as they were proposed in [ |J|||l4|] - yield the lowest 
energy expectation values. Here, a strong coupling ap- 
proach (a" 1 — > 0) to take into account singlet-triplet 
interactions is applicable. With this method Momoi and 
Totsukaliil found evidence for plateaus in the magnetiza- 
tion curve at M = 1/4 and M = 1/6. However, they did 
not find plateaus at smaller magnetizations (M = 1/8 
and M = 1/16), "since the mechanism to stabilize these 
plateaus is not yet clear" - as they say. 

We think that this failure has a simple explana- 
tion: The singlet-triplet configurations on the Shastry- 
Sutherland lattice (cf. e.g. Fig. |(c),(d) for M = 1/8) 
are unfavorable, since the formation of tripl ets o n next- 
nearest neighbor sites costs energy [cf. e.g. (|2.9|)1. Con- 
figurations - like K 2 in Fig. |](b) for M = 1/8 - with 
well separated monomers (unpaired spin-up states) yield 
a lower energy as long as a is not too large (a < a c {M = 
1/8) = 2.3). 

If the coupling a - realized in the compound 
SrCu2{BOz)2 - is indeed below this value, the exper- 
imentally observed plateaus at M = 1/8 and M = 1/16 
cannot be associated with singlet-triplet configurations 
on the Shastry-Sutherland lattice. 

In order to find the correct magnetic order at low 
magnetizations M < v/2N and a < 2.3 a more gen- 
eral ansatz for the ground state configurations is needed. 
This can be constructed in terms of monomers at fixed 
sites Xi . . . x„. The spins on the remaining sites form 
an antiferromagnetic cluster, the ground state energy of 
which depends on the fixed positions of the v monomers. 
Therefore, a specific distribution of monomers charac- 
terizing the magnetic order is given by a minimum of 
the ground state energy of the corresponding antiferro- 
magnetic cluster (cf. e.g. Fig. |](a),(b) for M = 1/8). 
We expect that for small values of M - in particular 
in the sectors with a finite number v of monomers, i.e. 
M = v/2N — > for N — > oo - the singlet- triplet config- 
urations on the Shastry-Sutherland lattice are dominant 
again (for a > a c (M = 0) = 1.4). Each of the (N - v)/2 
singlets lowers the energy by — 3a/4 whereas each of the 
few (v/2) triplets costs energy a/4. 

It should also be noted that the frozen monomer ap- 
proximation becomes better and better for M — > 0, since 
the antiferromagnetic clusters cover more and more of 
the whole lattice. 



Finally, we have also studied the formation of plateaus 
in the magnetization curve of the Shastry-Sutherland 
model. 

We looked for the a-dependence of the lower and upper 
critical fields as they follow from exact diagonalizations 
on finite clusters with N — 4 x 4 = 16 and N = 4 x 6 = 24 
sites. All the critical fields are almost a-independent for 
a < 1.2, but change rapidly above this value. Indications 
for the opening of a plateau are visible for M = il/ir4 / 6 
supporting previous results with other methods .EjEj 

The situation for M = 1/8 appears to be more subtle. 
The lower critical field has a pronounced minimum at 
a = 1.5. Here, the finite-size effects arc rather small. In 
contrast the upper critical field reveals a strong finite-size 
dependence. Computations on larger systems are needed 
for a reliable estimate of the thermodynamic limit. 
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